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ABSTRACT 


A quasi-geostrophic omega equation for large-scale vertical veloci- 
ties in the atmosphere is solved at the earth's surface, 850 mb, 700 mb, 
500 mb, and 300 mb for a grid of 1977 points per level covering the major 
portion of the Northern Hemisphere. The effects of terrain and surface 
frictional stress on the vertical velocity at the lower boundary are 
investigated in some detail and applied at the earth's surface. The 
three-dimensional integration of the resulting equations excludes areas 
"inside" the terrain. The mathematical model is programmed for the 
Control Data Corporation 1604 digital computer to utilize the operational 
output of machine-processed data and analyses produced by the Fleet 
Numerical Weather Facility (FNWF), Monterey, California. The program 
allows a wide variety of modes of operation and yields a solution for W 
in five minutes. 

An unusual weather situation in January, 1962 is investigated and 
the results, though tentative, show that the correlation between computed 
vertical velocities and the observed weather is excellent. 
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butions and assistance in this investigation. 
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he. Introduction 

The value of vertical velocities in the atmosphere for understanding 
the basic weather phenomena associated with the large-scale features is 
gradually becoming more fully appreciated by synoptic forecasters [1]. 

The purpose of this investigation is to provide a means of computing these 
vertical velocity fields from routine weather charts and to examine the 
Structure of these fields. Among the various methods for such determina- 
tions, the omega equation, derivable from the vorticity equation and the 
first law of thermodynamics, appears to afford the best means of computing 
large-scale vertical velocities from synoptic meteorological observations. 
Thus, the computer program developed for this omega equation is adaptable 
to operational use in a numerical forecasting center as well as for basic 
research. 

The data used in this program are produced by the Fleet Numerical 
Weather Facility (FNWF), Monterey, California, essentially as follows [2]: 

(a) Northern Hemisphere weather data are processed automatically by 
digital computer to yield reliable listings of surface and upper-air data. 

(b) A Northern Hemisphere surface map is analyzed by digital computer 
and 1000-mb height and temperature charts are produced. 

(c) The 850-mb, 700-mb, 500-mb and 300-mb height and temperature 
charts are then constructed with hydrostatic consistency by building 
upward from the 1000-mb level. Winds and off-time reports are also 
included. 

(d) These surface and upper-air charts cover the entire Northern 
Hemisphere to the equator, and are produced at 0000 GCT and 1200 GCT 
daily. Further details are available at the Fleet Numerical Weather 


Facility. 


7p Method of Solution 

On the basis of the available data described above, it was decided 
to arrange the computational form of the omega equation so that only the 
data at a given constant-pressure level would be needed to calculate the 
parameters required in the iterative Scheme for that level. This necessi- 
tated the use of unequally spaced pressure intervals in computing 
vertical derivatives, but gave the model some flexibility if additional 
levels became available or were changed. 

The omega equation is derived from the vorticity equation (1), and 
the thermal equation (2) for adiabatic conditions by eliminating the 


local time derivatives jointly between the two equations. 


oz oF _ dw du _ dw du , 
Seen) ap = 1S +e dx dp ~ 


£2 4VeVT-wT =O 3; T=-+—, (2) 


where the various symbols are defined in Table of Symbols, page vy. 
The geostrophic approximation for vorticity has been made in the first 
term of equation (1). Operating on (1) with - and on (2) with 


-7 } and adding results in a form of the omega equation 


(3) 


ee me 


3 dP ox OP 
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The terms involving the vertical advection of vorticity and the 
twisting terms are usually of opposite sign and, moreover, small in com 
parison with the others remaining, and will be neglected as a first 


approximation. Of the terms which arise from the horizontal derivatives 


ang a, 
9 f 2. De . 

V (twleTV w+ wT +20T°Vw . (4) 
the last two on the right are also considered to be small in comparison 
with the first and will be neglected. The neglected terms may be 
incorporated along the lines used in Stuart [3], to afford a solution 
to the more complete equation (3). It is to be noted that 9 will be 
allowed to vary in the horizontal as well as with pressure in this scheme. 
To evaluate partially the effect of the variable &@, the computer program 
has been designed to utilize alternately horizontal and vertical means of 
the static stability as well. 

Neglecting the last two terms in (4) and substituting in (3) gives 
the reduced omega equation (5) 


o vin Pla 2% _ yyevt)- PL 2(vev7) =0 5) 


R Sp 
where, in the second term, fr has been replaced by i to meet certain 
consistency requirements as described in [4] and [5]. Again, to evaluate 

the effect of this term, the computer program has both versions Ty and t, 
incorporated. 


The third and fourth terms of (5) can be expanded using the geo- 


strophic vorticity ; b= <i var: : 


the hydrostatic equation, and the equation of state: 


iv (V- > - - mia v7) = 


ss 2. fe = Zz 
= |-V 5(2,7)-E5(%0) +J(Z,V 7) 
Substituting this expression into (5) and rearranging gives the computa- 


tional form of the reduced omega equation: 


Vor Pile = = Slvr aT) s (2 VT)- 17,0) (6) 


mv) 


The basic parameters involved above are Zz and TT - At a given 
level, the right side of (6) can be calculated from data at that level 
only. The last term on the right in (6) is generally considered to be 
of major importance and is used in the approximation made by Wiin-Nielsen 
[6]. 

The grid used is the Fleet Numerical Weather Facility Northern Hemi- 
sphere octagon of 1977 points centered on the pole and extending to about 
latitude 10N. 

The omega equation (6) is to be solved at the four levels of 850 mb, 
700 mb, 500 mb and 300 mb, with the lower boundary being the surface of 
the earth. The upper boundary is 100 mb where it is assumed that W=O , 
and the lateral boundaries are near 10N where it is also assumed that 
pn=O© . 

To solve (6) numerically, the following finite-difference approxima- 


tions, based on the adjoining five-point array, were used. 
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I ® 
O8 
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A.=As-A, 5 AYAZA-Ag y BpA=|(A,AY Wn SF, 
A A=A,+A-2Ao 5 AA =A,+Ag-2A, ; _ 
VA=MA+MA » VA=E7A 5 

‘i (A,B) = (B.A): (4,8) - (A,A):(.,8) » J(A,B)= m5 7 (A,B) - 
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The portion of the term involving P Vm 1S approximated [4] by a finite- 
C 


difference analog for unequal intervals 
, 2 
: OW = aie ZA ww (8) 
Z 
OP Ayp Ap 


where Ae ee AxP LO + AvP Wee 
Dip Ay, p 


and AP ine and Ap are defined in figure l. 
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The static stability 


a eS (9) 


is equal to the difference between the dry adiabatic lapse rate 
bp = RT/ pcp in pressure coordinates and the observed lapse 
rate Yp = OT/ 2 P . This form of the static stability parameter 
corresponds to that suggested by Gates [7] as "one of the more useful in 
modelling the large-scale disturbances in a troposphere of constant 
s@apility . 

Since it is evident that the value of 0 may approach zero or be nega- 
tive, the magnitude of @ was controlled so that the second term of (6) 


remained positive and within certain bounds. This was done by requiring 


that the lapse rate not exceed 7/8 of the dry adiabatic lapse rate, namely, 
an, NYE 
{p — 3 (op or 
* 7 ~3 
gg => (P,p i ; (9A) 
pp 2 
Comparison of the true O with PD z for a few Synoptic cases 
has shown that this limiting criterion is very rarely, if ever, required 


when the input data are vertically consistent. 


The second term of is approximated by the finite difference analog [8] 


aT) = Dole =| Pe MA Te + SEA Ts 


OP Ap MP AP AYP. vl (10) 


where A, and Ne are defined in figure l. 

Sigma is expressed in units of degrees Kelvin per mb. [In order to 
calculate sigma at 300 mb, temperature values at 100 mb are required. 
These data are not presently available from the Fleet Numerical Weather 
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Facility on a routine basis, and hence temperature values at 100 mb were 
approximated by either of the following schemes: 

(a) The temperature at 100 mb was set equal to the 300-mb temperature. 
This approximates mid-latitude conditions and implies that a "tropopause" 
ex1sts Somewhere between the 300-mb and 100-mb levels, the specific press- 
ure depending on the fitted parabola, 

(b) A mean 100 mb temperature field is used based on a climatic 
summary. Since the synoptic case chosen for study was in January 1962, 
5-year mean 100-mb temperatures for January (1949-1953) were used. 

It is realized that neither of the above schemes is entirely satis- 
factory; however, the computer program iS so arranged that synoptic analy- 
ses of 100-mb temperatures can be used when available on an operational 
basis. 

In order to ut1lize established numerical sub-routines available at 
the FNWF as well as to minimize computational time, fixed-point fractional 
arithmetic was used excluSively in the computer program. Moreover, all 
parameters had to be scaled to have magnitudes less than one. The follow- 


ing standard FNWF scaling conventions were uSed: 


I 


eZ (mb sec!) n=8 at present 
ae | a 
=£-2" (see") 


se) na p-2" Gee! 
=p:2" (em) Z=7 
g 


2°? (am) (11) 
Z 


3 Uo we 
| 


= Pree (mb) = G2? (°K mb) 
= re? + ee (°K sr *C) 


Subs tustutdmgen(7), GBdsy C9d. (10) swamd (11) iitiomG6) cives the 


Oo 


scaled difference equation that was programmed for the digital computer: 


where B. = B..2 = “ee py zl ; Px AL | ° 


(> 
lll 


meeting a Far st, Tenmeonen ie micaima@E 2) a, 


Second term on right in (12) 


equation (12) becomes 
Vodx — Bet = Ce -Ca HECK K=L2 6? LOD 
K K KK J SE re Ii ) ° 
These equations were solved by an accelerated Liebmann relaxation 


; yo 2 ‘ A 
iterative scheme where the (UW+/ ) iterate at the point Ce, 4 ) is given 


by 


Sar = ee (14) 
4+B 
Here the residual R at each Step is: 
ua ZBAav WoL 
K aa VY UW, - Bi —Cr + Cyr ‘ (15) 


The four pressure levels were successively scanned beginning at 


850 mb and ending with the 300-mb level, this constituting one SCAN count 


9 


of 4(1977 )= 7908 points in the 3-dimensional mesh. 

Considering (13) as a system of Helmholtz equations, certain condi- 
tions regarding the coefficients must be fulfilled for convergence 
('Arnason [9]) as follows. 

(a) Be > O (16) 

(b) The sum of the coefficients of W_x-, and Wy, In CA im. (13) 


must be less than or equal to one, ineluding round-off errors 


' AN D : Not | ae . . 
Bi Ax? i Ae oe (17) 


The first requirement (16) is easily met since Vk PE, (see (9A)). The 
second requirement iS sensitive because of the great range of the co- 
efficient Bi ,» Which can vary from a very small positive value (arbitrar- 
ily set ee to values of the order of 2000. Both of these require- 
ments are incorporated into the program. 


The convergence criterion € is defined as 


for all points of the mesh. The € is chosen so that no meteorologically 
significant changes will occur in the omegas if iteration were continued. 
The smallest value of interest was taken as lone ne eee which is approxi- 
mately equal to 0.1 cm ose at 650 mb: 


A number of different values of € , ranging from oer to (1/6400) 


ee nb eae » were tried to determine an optimum © . The final opera- 
me = 
tional value adopted was Gre = — Ne mb sec i 


With this criterion, Cos , convergence is reached, with the synoptic 
cases studied, in 14 to 15 SCANS as defined previously. For comparison 


purposes, convergence for Cise= a. oc mb om is 
oO 


‘ Loe 


obtained in about 26 SCANS, 


The over-relaxation coefficient IN used was determined experi- 
mentally by finding that A for which convergence was reached in the 
least number of SCANS for a given € . This optimum is very 
sharply defined and was found to be 1.414. A A differing from this 
value by roughly 10% increased the SCAN count by about 25%, and a change 


in value of only 0.1% changed the SCAN count by one. 


De The Lower Boundary Condition 

The lower boundary condition was of special interest in devising 
this model [10]. In most previous studies, the given or calculated lower 
boundary values of omega were arbitrarily assigned to the lowest constant- 
pressure level in the model, usually 1000 mb, and then the family of 
equations was integrated with no regard to the presence or absence of 
mountainous features. The intersection of any pressure surface with the 
terrain constitutes an internal boundary within the mesh and is treated 
as such in this model. It is noted that, in order to find the inter- 
section of the terrain and a given pressure level, the pressure at the 
terrain Pr must be known. 

The input data at the various pressure levels (K) consist of Dy 


values and temperature Ty in degrees C, The definition of D is: 


Do 78 fie . 


The D-+- value at the height of the terrain Zy is interpolated from the 


input data Dy by the principles of altimetry ee 


ei 
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AyD = Dy -D, 
A, Zp = Zp.-ZLp, (18) 
My.Zp = Zp. fa Tae 


A.D Dy = 1D er We = y= AyD 2 


Vv 


\\j 


To find ie D , the specific temperature anomaly Syy =A D/A Zp 
memcalculated. This 5j, value is the mean for the column Pv to Pi 
Applying it to the portion from Pu to er 

nN _ a 
Substituting in (18) gives 


Dp. = Du=Suri( Zp. - Zr) . (19) 


} A ae Sf 


Having interpolated D~- , the pressure altitude is then 
Zp; = a Die =: 
and i can be obtained from 2P; Since there iS a one to-one corre- 
spondence of values. 
For the last step, an interpolating polynomial is developed that 
gives the terrain pressure with accuracy within one millibar. 
The polynomial is 
7a : wee 
aa oS 16 *’). (c- , 20 
P+ y(at by") +U(c+du*) , (20) 


ls: 


where a = 724.64 c = 597.37 
b = 3.5562 d = 3.1272 
y = |-u 
u We. 272 


L5-/07 


The scaled polynomial is then 


Pr =2°|y (44 697) +0(@4+d0")] , 


where a = 4.4229,1075 

Pee) 2-47 22 aloes 
C = 3.6461.10 5 
d 3.0539.10 

A ce 
u =uUu'Z 

roe = 

=y-Z 
ee . 
Ge e2 pr = 0.25 7a 
0.5 7220 

Y= ena ey oe 


The pressure at the terrain is calculated on a synoptic basis and is used 
in devising a marker field which controls the iterative scheme so that | 
all points of the mesh which "lie" within the terrain are not relaxed 

and constitute an internal boundary. At the surface of the terrain a 

value of W,.,. is computed as described later. This W,,,value is assigned 

to the same grid point at all greater pressure levels, i.e., those within 

the terrain. This substitution is illustrated schematically in figure 3. 

Note that under this scheme, the Laplacian of W at (c) involves the 

uses of (Oro dir] , which generally more closely approximates the 


value of Wy than 1f zero were used. 
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4, Omega at the Lower Boundary 

The pressure at the terrain Pr is also used to calculate the omega 
at the lower boundary, W,.,. The combined effect of the frictionally- 
induced vertical velocity and the vertical velocity due to the slope of 
terrain will be calculated for the volume of air bounded by the surface 
of the earth at pressure ane and at the top by a Pages e me where the 
surface frictional stress is assumed to vanish. The net inflow or con- 
vergence into this volume results in a vertical velocity Wy,through the 


upper boundary. If the equation of continuity 


7 ae = AAV (21) 
Op 


is integrated through the friction layer, Pv — Pro , utilizing the 


wind V , 


vaH[ge roi ]e+ {fee 28]j. om 


= Bi (reop) + gf 8 he], 


where i is surface stress at Pe 

and ie a ty mC ail Pie 

The first term on right in (23) is usually small in comparison with the 
other terms and is neglected. Setting Wy, -Wy = Wpe meehe 


frictionally induced vertical velocity (23) may be written as 


A) fs 


a $) ates _ ci | . 
7 Oy Ox 


(24) 


If it is assumed (Lettau [12]) that the surface stress is opposite 
the direction of the wind, V- , making an angle & with the isobars, 


(see figure 4), then 
° sna AL one 2 Ay oom AL > 
a coker COS ) loy = lo Sing o 


The surface stress is further assumed to be proportional to the square 


of the wind at the top of friction layer in the manner of Lettau [11] 


ee 
and Cressman [10], i.e., ie, = Cp Fhe VAS 
and Tex = Cp Pic U..Vio : Tey =Cp Lo Vio Vic 7 (25) 
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Then from (24) 


_ 36,6, CoWMec)  3(CovieVe) 
a 7 Oy Ox | 


where the local variation of density has been neglected. Expanding this 


equation 


=a | oe cae yy) \ a 
WD) = = g GI CVs (cies ee vs) a U is Een) ~ hr (, CrVic) 
= 7 


rf 
| 
Oy ge 


°( 26) 


The last two terms on the right in (26) were assumed small in comparison 
with the first and were neglected. The surface friction omega can then 


be written as 


We = - PRE COV ie a 2 (27) 


where the density has been replaced through the equation of state. 
Next, the vertical velocities due to varying terrain height will be 
considered. The value of omega due to the terrain Slope@--may be written 


in the form 
dus) za 
()- =(s2 = \V/ : ra) an > 
PS Kdper Vocly © “te VPs — 
For a large 3-hourly tendency 
(ep = 5 mb per 3 hr 
ee 7 | 
and will be neglected. Then 
a =V_o°Vpr ° oe 


i 


and omega at the lower boundary is 


Ver -F ae pVic Fie : (30) 


The geostrophic wind at the lower boundary Vio is required in (30) in both 


terms. The geostrophic vorticityh,o is also required and is based on Vio 
V.o= kx ZV D)._ (31) 
Lo £ ie 


Here it is noted that Vp_D must be determined at terrain pressure. Thus 
it is necessary to ''reduce" or extrapolate the D=- values at points sur- 
rounding the central point. These extrapolated values will be called 

(Dr), . The extrapolation will be done on altimetry principles using 
the previously calculated Sy, values for the layer in which the terrain 


lies. For example, in figure 5, at point ( érL ), the extrapolated value 


is 


(Dr)p, etl (D-) iri * (Sux etl 


Zp.- f pr, det | (32) 


It follows that the gradient of D in the x-direction may be written 


(De)p int ~ (Dr) p é-1 = ID, Lel Dr, i-1| a Z od] (Sudie-(Sudéel 


(33) 
+ (Sudi-nZpr,é-1 ee (Some eee Ney P 


Similar considerations hold for the y-component. Both equations (32) and 
(33) are used in making calculations of Vig and eieor use in equation 
(30). 

The T+ in equation (30) is also calculated by altimetry principles 


[9]. By definition, the specific temperature anomaly is 
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S = tae, (34) 
‘Tp 


where ea the pressure temperature, is the standard-atmosphere tempera- 
ture at a given pressure altitude and T is the temperature at a given 


point. Solving (34) for T and applying to this model gives 


The pressure temperature There given by 
za 6 
Py = Wess Cle f Pr (36) 
where it = 15C is the temperature at Z = O 
Po Pp 
F ee 
and a = ©.5 10 °C per cm. 


Combining (35) and (36) yields 
is = (TR. —-OrF pz)({-Ser) : (37) 


Substituting (31), (32), (36) and the geostrophic vorticity into (30) gives 


the computational forms for W+ and O) = 
9 Tip). 027 
Pro (Dr)o5 Pr] > oe 


a: —*Cp:|Kx WV, D+] VACD) 
FPR Giles = Q' Zpr)(1t Su.) 


(32) 


Using the finite difference analogs (7) and the standard scaling conven- 


tions (ll) gives 


Sy = Fe Pl(b)p Pr] “0 
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We=- = - a (41) 
a £ sind: d?-R (Tp. -& Zpr)(it Su) 
A AY a A 
The , ,is then (Oi = Ud 7 +e > 


The computational time to obtain (3, is about 40 seconds. 

The coefficient of drag Cp in (30) is the same as that defined by 
Cressman [13], and values of Cp were obtained from the Joint Numerical 
Prediction Unit, United States Weather Bureau, Suitland, Maryland. 

The height of the terrain Z7 is an area-smoothed value, and for the 
Synoptic cases presented the heights used were obtained from the FNWF. 
The terrain heights have certain deficiencies due to rather heavy smooth- 
ing. They are being replaced with values averaged over a l-degree 


latitude/longitude area, as determined by Berkofsky and Bertoni [14]. 


oe Some Features of the Computer Program 

The Control Data Corporation 1604 digital computer was used in this 
investigation. It has a core memory of 32,768 words of 48 bits each, 
with a typical add time of 7.2 microseconds. 

The initial "guess" field for omega at the various pressure levels 
is obtained from the W,,, value previously computed. It was expected that 
the influence of a particular (Uat the lower boundary would decrease 
nearly exponentically with height above the terrain, hence the initial 

(c) Pre: 
guess field was computed by (Ol ee eae ) 
where t = 1 for the first pressure level above the terrain, t = 2 for the 
second, etc. This procedure accelerates the convergence of the relaxa- 
tion scheme appreciably, as compared to an initial guess field of zero. 


Provision is made in the computer program to determine the influence 
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of the lower boundary on the omega fields in various ways. 

(a) The terrain-induced vertical velocity W- and the frictionally- 
induced velocityWemay be calculated separately or jointly. 

(b) &,. may be set equal to zero. Here the terrain is still present 
and the lower boundary follows the variable height of the surface of the 
earth. 

(c) The computation scheme can be altered to test the effect of 
inserting the W,;,values at an artificial lower boundary, Calculated on 
the assumption that the height of the terrain is zero. The pressure at 
this lower boundary is not constant and corresponds to the sea-level 
pressure map. This Simulates closely the approximations normally used 
[3], [10], [13], wherein the lower boundary values are arbitrarily 
assigned to the lowest constant-pressure level of the model. Thus, they 
are not only applied at an incorrect altitude, but ignore the presence of 
internal boundaries due to the intersection of pressure surfaces and 
terrain. 

(d) The forcing function Cz of equation (3); whieh anvolves the rl- 
and T-fields responsible for the large-scale vertical velocities, can be 
set equal to zero and a Solution obtained with only the lower boundary 
values Wi, present. This essentially shows the decay of the influence of 
the Gd), -, with height, subject here only to the static stability variations 
in the atmosphere (See equation (6)). 

Examples of the various possibilities (a) to (d) will be given in 
the next section. 

Many modes of operation of the computer program are possible and 
allow a wide variety of side products to be obtained if desired. The 
more important of these are listed below. 
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D, at 5 levels 
T. at 6 levels 


Gg eatea levels 


Wio 

f at 4 levels 

Ww at 4 levels 
All of these can be obtained in computer-analyzed form on machine-drawn 
and labelled charts. All the charts of this investigation were produced 
in this manner except cross sections. 

It is to be noted that on the charts of Wo wherever the terrain 
projects through a pressure level, the value +1lll is printed at that grid 
point, and delineates the internal boundary. 

The computer program has "smoothers" incorporated to control irregu 


larities in higher-order finite differences of the input data. These 


smoothers are of the form 
ea Z 
ee eet KW ALS 


Na 
where A is the smoothed data, 
A is the unsmooth data, 
and KK is the smoothing coefficient and is equal to 1/8. 


2 Ke jan Aoa/s > 
This smoother is applied once to the Y De 5 V IK 5 and V y OF Te) 


fields. 
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Oe Discussion of Results 

The static stability 1s a bulk value, since it is derived by fitting 
a parabola to temperatures from three pressure levels of hydrostatically- 
checked analyses. The computed Gx vary as expected with regard to the 
relative stability characteristics of various air masses at lower levels 
and with respect to the proximity of the stratosphere at upper levels. 
The observed rapid variation of the static stability may not justify the 
neglect of the derivatives of Gx in certain areas. These effects are to 
be included in future computations. Some typical charts are displayed in 
Appendix C, pages 136, 139, 142, and 145. 

The influence of the lower boundary on the vertical velocity field 
was investigated, and some typical charts of the results are displayed in 
figures 6 to 10. Note that the decay of velocity with distance from the 
surface of the earth is very rapid, and even large values of w,. (at 
levels below 850 mb) do not exert a significant influence beyond the 
500-mb level. 

In many models, the vertical velocity at the lower boundary 1s 
artificially assigned to the lowest pressure level of the model. The 
present model was constrained to operate in this mode for a few tests, 
one of which is presented in figures 11 to 15. Inspection shows that, as 
expected, there is little disagreement in areas where the terrain is near 
zero height, but large discrepancies arise at those points in the vicinity 
of prominent terrain features. Areas where this disagreement is the 
largest are obviously the Cordilleras of North America and the Himalayan 
Plateau of Asia, The transparent overlay at the end of this report may 


be uSed on figures 6 through 10. 
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Figure 6. 
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ie Mereical-Velocity Computations for a Particular Synoptic Situation 

The California snow storm of January, 1962 was selected to provide 
data for vertical-motion analyses. It was considered that this storm 
might provide results which would show the development of a situation 
which was inadequately forecast. An extenSive synoptic investigation of 
this storm has been made by Ikegami [15]. 

The required input data of D values and temperatures at 1000, 850, 
750, 500 and 300 mb were obtained from FNWF, Times chosen for analyses 
were 1200Z 17 January, 1200Z 20 January, O000Z 21 January, 1200Z 21 Janu- 
ary, 0000Z 22 January, 1200Z 22 January, 1200Z 23 January, and 1200Z 
26 January. 

Appendix A contains extracted portions of hemispheric prints which 
measure 30 x 26 inches. The extracts are the same scale as the originals. 
Since various cross sections are provided, a line delineating the cross 
section path is indicated. Also, grid points three and four are under- 
lined. The Monterey Bay area lies between these two points, as shown by 
the continental outlines drawn on the chart on page 44 of Appendix A. 

This extract gives the pressure at terrain height for Q0O00Z 22 January. 
If it is desired to use the transparent overlay on these charts, shift 
the overlay to the left one grid point for proper display. 

As indicated earlier, wherever terrain projected through any pressure 
level, the value +111 was printed at that grid point. For the maps pro- 
vided, only the 850-mb level is affected. Since a point within the 
terrain was not relaxed in computing the omega fields, the terrain value 
of omega was substituted at that grid location. 

The printed values of omega multiplied by 1077 give units of mb/sec. 
Vorticity values multiplied by i0~° give units of peoia. 
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a. Analysis for 1/7 January 

(See Appendix A, pages 45 to 54 ; Appendix B, page 106) 

In general, California was under the influence of subsidence associ 
ated with a high located to the west. However, in the northwestern United 
States a wave developed on a Stationary front oriented from northwest to 
southeast. This wave continued to develop and the associated cold front 
brought heavy rains into California. 

Rain began in the Monterey area about 1500Z 19 January and continued 
through 1300Z 20 January. Frontal passage occurred about 1000Z 20 Janu- 
ary. The pressure rise was small and no appreciable increase in pressure 
was noted for 48 hours, due to a surface trough which remained in the 
area. 

Program outputs of the region where this development occurred depict 
ascending air at all levels. Further, relative vorticity analyses show a 
maximum located northwest of the centers of rising air. 

Appendix C, page 119, shows the presence of a blocking high at the 
500-mb level. The block's orientation and the downstream cold low became 
the controlling factors in the movement of the surface wave. 

b. Analysis for 20 January 

(See Appendix A, pages 54to 59; Appendix B, page 107 .) 

By 1200Z the front had passed Monterey, but intermittent showers 
continued to about 1800Z,. The center of ascending air was located north- 
east of the Monterey area with higher vertical velocities at the upper 
levels. Strong subsidence to the northwest was apparently moving east- 
ward and would soon dominate the local weather. 

But, as seen in Appendix C, page 121, the previously-mentioned block 
was now fully developed. The downstream low became quasi-stationary and 
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the area of subsidence remained to the west. 
c. Analysis for 21 January 

(See Appendix A, pages 59 to 77; Appendix B, pages 108 to 110.) 

The 500-mb chart for 0000Z (Appendix C, page 128) shows essentially 
the same situation as existed 12 hours earlier. The temperature of the 
cold low was near -40C. Moreover, this tongue of cold air extended down 
to the 1000-mb level. However, the ascending air was now located south- 
east of Monterey, and higher vertical velocities were found at lower 
levels. 

The 1200Z analyses show that the center of maximum relative vorticity 
and ascending air were coincident at low levels, i.e., through 850 mb. 
Using the terrain boundary on the 850-mb chart, it appears that topography 
restrained the lower portion of the rising column, shearing it from the 
main cell. 

Rain commenced locally at 0600Z 21 January. Snow pellets began 
falling about 1400Z and snow showers began at 1615Z. Throughout the 
period ending 1000Z 22 January, snow showers, snow pellets and rain were 
frequent. 

d. Analysis for 22 January 

(See Appendix A, pages 77 to 95; Appendix B, pages 110 to 115.) 

The significant aspect at QO00Z was the change in direction of move- 
ment of the cell of rising air to the southwest. Vertical velocities 
decreased near points three and four at the 850-mb surface. Appendix C, 
pages 136, 139, and 142, shows the extent of low-level instability. How- 
ever, to the west, subsidence was clearly increasing at all levels and 
now appeared to be moving to the east. 

The 1200Z results show further movement and intensification of 
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subsidence as it proceeded into northern California. The ascending air  & 
continued moving southwestward. The surface low, after having reached 

a position in southeastern Nevada, moved southwest and weakened west of 
Lower California. At 1200Z 22 January, it was situated over the Los 
Angeles area. 

eo Analysis for 23 January 

(See Appendix A, pages 95 tol00; Appendix B, page 115) 

These 1200Z charts, 24 hours later than the previous analyses, show 
points three and four were under the direct influence of descending 
motion, save for a small area of rising air due to terrain effects near 
point four. The velocity of the ascending air to the southeast generally 
decreased and the cell appeared to be weakening. 

In Appendix C, pages 146 to 156, the shift of cold temperatures and 
the trough towards the southwest and the weakening of the block at 500 mb 
are depicted. 
ae Analysis for 26 January 

(See Appendix A, pages 100 tol05; Appendix B, page 116.) 

The entire West Coast was under the influence of subsidence and 
little evidence remained of the column of rising air. Appendix C, 
page 157, shows the trough had resumed normal west-to-east movement. 

2. Conclusions 

The preceding analyses suggest that had this vertical velocity data 
been available, the forecast of fair and sunny conditions for 21 January 
would not have been issued. The vertical motion fields, together with 
temperature analyses, indicate the impending storm. 

The omega fields were compared to the surface analyses as trans- 
mitted by facsimile, and indicate that the direction of vertical motion 
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might have been inferred from conventional map analyses alone. However, 
the magnitude of the velocities of ascending air in the surface trough 
could hardly have been appreciated from such maps. 

Lack of time precluded computing vertical velocities for other 
synoptic situations. Also, as interesting and encouraging as these 
indications are, it should further be pointed out that the results dis- 
played are for an area of abundant data. It remains to be seen how uSe- 
ful and accurate vertical-velocity computations are in areas of sparse 
upper-air data. 

The outputs of this program can form the basis for other investiga- 
tions in the field of numerical forecasting. Presently, the vertical 
velocities obtained here are being used to compute the transformation of 
potential energy to kinetic energy on a hemispheric scale. 

With prognostic D and temperature fields, prognostic vertical 
motions can be computed. The addition of hemispheric dew-point analyses 
and prediction would permit forecasts of cloudiness and precipitation. 
The contribution of the release of latent heat to the stability parameter 
could then be included. The obvious application of the omega fields to a 
multi-level forecast model is another possibility. 

On the hemispheric charts, numerical shading between plus and minus 
areas of vertical motion imparted an appearance of cloud patterns. 

Since facsimiles of TIROS IV cloud pictures were available, two samples 
were chosen for comparison. 

Appendix D contains the TIROS pictures and the applicable omega 
fields. Page 159is the 0501Z 17 May 1962 TIROS observation; an apparent 
vortex is indicated near 35N 163E. Here the vertical motion charts lag 


the observation by seven hours. 
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Page 165 is a 2220Z 29 May 1962 picture. A vortex is located near 
42N 100W. The omega charts are for a period about two hours later. 

The observations and representations are well correlated. Perhaps 
further study would enable the cloud type to be correlated with the 


vertical-velocity field with respect to sign, magnitude, and level. 
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APPENDIX A 


EXTRACTS FROM HEMISPHERIC PRINTS OF 


OMEGA AND RELATIVE VORTICITY FIELDS 
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APPENDIX B 


CROSS SECTIONS OF OMEGA FIELDS 
PLUS TEMPERATURE, RELATIVE VORTICITY AND STABILITY 


FOR O0000Z 22 JANUARY 1962 
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APPENDIX C 


HEMISPHERIC CHARTS OF D VALUES 


AND TEMPERATURES PLUS STABILITY 


FOR 0000Z 22 JANUARY 1962 
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